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Numerical solution of damped forced 
oscillator problem using Haar wavelets 
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Abstract 


We present here the numerical solution of damped forced oscillator prob- 
lem using Haar wavelet and compare the numerical results obtained with 
some well-known numerical methods such as Runge-Kutta fourth order clas- 
sical and Taylor Series methods. Numerical results show that the present 
Haar wavelet method gives more accurate approximations than above said 
numerical methods. 
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1 Introduction 


During the last few decades considerable efforts have been made using 
wavelet, towards the development of computational methods to solve numer- 
ically linear differential equations encountered in various fields of science and 
engineering. Wavelet analysis is a new branch of applied science. Wavelet 
methods are applied to find the numerical solution of problems related to 
science and engineering. In the last recent years, wavelet methods have been 
attracted the great interest of researchers of physical and mathematical sci- 
ences and many research papers were published in these fields. Recently, 
many researchers have used Haar and Daubechies wavelets to find the nu- 
merical solution of differential and integral equations. Before, the discovery 
of Haar wavelet, Daubechies wavelets were used in many published research 
papers for numerical solution of differential and integral equations. 

In 1910, Alfred Haar [4] discovered a new wavelet known as Haar wavelet. 
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Among all wavelet families, Haar wavelet is most simple, accurate and effi- 
cient. It attracted, the interest of many researcher in the field of engineering 
and science. Haar wavelet has been used in wide variety of numerical meth- 
ods developed for numerical solutions of differential and integral equations. 
Here, we present a survey of such methods for differential and integral equa- 
tions. Chen and Hsiao [3] applied Haar wavelet method for solving lumped 
and distributed-parameter systems. Hsiao [6] used wavelet approach to time- 
varying functional differential equations. Razzaghi and Ordokhani [15] used 
Haar functions for variational problems. Ohkita and Kobayashi [13] applied 
rationalized Haar functions to solve linear differential equations. Cattani [2] 
suggested use of Haar wavelet splines for numerical solution of differential 
equations. Lepik [8, 9, 10, 11, 12] used Haar wavelets for solving differ- 
ential and integral equations. Sunmonu [18] presented wavelet solution for 
second order differential equations with maple. Hariharan and Kannan [5] 
presented an overview of Haar wavelet method for solving differential and 
integral equations. Kouchi et al. [7] presented numerical solution of homoge- 
neous and inhomogeneous harmonic differential equation with Haar wavelet. 
n [16], Quasilinearization technique and Haar wavelet operational matrix 
method both are applied to find the numerical solution of fractional order 
nonlinear oscillation equations. Also, Solutions of fractional order force-free 
and forced Duffing-Van der Pol oscillator and higher order fractional Duffing 
equation on large intervals are presented in [16]. 

In Section 2, we discussed damped forced oscillator. Haar wavelet method 
is presented in Section 3. Function approximation is presented in Section 4. 
In Section 5, we present convergence analysis of Haar wavelet method. In 
Section 6, the solution by Haar wavelet method is presented. In Section 
7, Runge-Kutta method for second order differential equation is presented. 
Taylor-Series method is presented in Section 8. Comparison of numerical 
solutions is presented in Section 9 and in Section 10, conclusion is given. 


2 Damped forced oscillation 


Oscillation means repeated motion of a particle or a body, when displaced 
from its equilibrium position. The classifications of oscillating systems are 
presented in Thomsen [19] and in Bhat Rama and Dukkipati [14]. The mech- 
anism that results in dissipation of the energy of an oscillator is called damp- 
ing. In mechanical oscillator, the damping may be due to (1) Viscous drag 
(2) Friction and (3) Structure. An oscillator to which a continuous excitation 
is provided by some external agency is called forced oscillator. 

Suppose a mass M attached to the end of a spring of stiffness constant 
S. One end of the spring is attached to a rigid support. Let the driven force 
acting on the system be F(t). At any instant, the system will operate under 
the influence of the following forces: 
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(a) Restoring force, F1 = —Sx where z is the displacement of the mass from 
the equilibrium position, 

(b) Damping force, F2 = —rdx/dt, where r is damping constant, 

(c) Driven force, F3 = F(t). 

The negative sign in the first two expression implies that both the restoring 
as well as damping forces opposes the displacement. By Newton second law 
of motion, we have 


M =—Sa—r— + F(t). (1) 


In this paper, we take special choice F(t) = 2(1—sint), M = 2kg, S =1N/m, 
r = 0.3Ns/m and x(0) = x'(0) = 0 as initial conditions, see Simmons [17]. 
The exact solution of equation (1) by using classical method is: 


2 
a(t) = e~ 9-975 (Cy cos(0.703118t) + C2sin(0.703118t)) + 24 OO aes Tc cos(t). (2) 


109 109 
applying initial conditions, we have Cy = —73 and Cy = - ee 


3 Haar wavelet method 


The Haar functions are an orthogonal family of switched rectangular wave- 
forms where amplitudes can differ from one function to another. They are 
defined in the interval [0, 1]. 


tl, @S0< 2, 
0, otherwise. 


k+0.5 = 
a and y = ‘ 


where a = z. B= 
Integer m = 23, (j = 0,1, 2,3, 4,....... J) indicates the level of the wavelet. 

k = 0,1, 2,3,.....,.m—1 is the translation parameter.Maximal level of resolu- 
tion is J. The index 7 is calculated according the formula i = m+k+1.In the 
case of minimal values, m = 1, k = 0 we have 1 = 2. The maximal value of i 
is i= 2M. where M = 2’. It is assumed that the value i = 1, corresponding 
to the scaling function in [0, 1]. 


1 0<t<il 
hi(ft)=4 07 4 
(6) ‘ otherwise. (4) 


Let us define the collocation points t; = ats where | = 1,2,3,...,2M 
and discredits the Haar function h,(t). 
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In the collocation points, the fist four Haar functions can be expressed as 

follows: 

hy(t) = [1, 1,1, 1], he(t) = [1,1, -1, -1], h3(t) = [1, -1, 0, 0], ha(t) = [0,0, 1, -1]. 


We introduce the notation: 


1111 
H,(t) = [hu(t), halt), hal) hal” =] 4) Gg |- (5) 
00 1-1 


Here H4(t) is called Haar coefficient matrix. It is a square matrix of order 4. 
The operational matrix of integration P, which is a 2M square matrix, is 
defined by the relations: 


Pa(t) = | Rae (6) 


0 
ty 
Ping i lt) = Pi n(t)dt, (7) 
0 


where n = 1, 2,3,4.... 


These integrals can be evaluated using equation (3) and first four of them 
are given below:- 


t— a, tela, B), 
Pit) ne a od t, te[B, as (8) 
0, elsewhere, 
g(t— a), te[a., 8), 
P,,2(t) = int > s(y- 1), te[B,Y), (9) 
. at tely, 1), 
0, elsewhere, 
a(t =) tela, 8), 
1 3 
Pra(t) = 4 HEA gly— Otel), ae 
imz (t — 8), tel, 1), 
0, elsewhere, 
ga(t — a)4, tela, B), 
P, a(t) = a2 (t B)? Xi (y i)? + went te[B,Y), (11) 
Fai? (t 6)? ant , te[y, 1), 


0, elsewhere. 
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4 Function approximation 


Any square integrable function x(t) in the interval [0,1] can be expanded by 
a Haar series of infinite terms: 


x(t) = )~ ajhj(t), ie{O} UN (12) 
i=1 
where the Haar coefficients a; are determined as: 
1 
ne i; x(t)ho (t)dt (13) 
0 
pl 
An = 2 | x(t)h;(t)dt (14) 
0 


where i = 2/ + k, 7>0 and 0<k < 2), xe[0, 1] such that the following integral 
square error € is minimized: 


€ =| [x(t) — ys a;h;(t)]?dt (15) 


where m = 2? and je{O}UN. 

Usually the series expansion of (12) contains infinite terms for smooth 
x(t). if x(t) is piecewise constant by itself or may be approximated as piece- 
wise constant during each subinterval, then x(t) will be terminated at finite 
m terms. This means 


a(t) = 3 ajh;(t) = ang? hin) (16) 


where the coefficients a7 and the Haar function vector h(t) are defined 
as: 


T= | 


adm AQ, A], AQ, vereeeeeee ,am—1] 


and 


5 Convergence analysis of Haar wavelet method 


Consider a differentiable function x(t) with 


|z(t)| < Ko, (17) 
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such that 
|x'(t)| < Ko, (18) 


for all te(0,1). Where Ko > 0 is a positive constant. Haar wavelet approxi- 
mation for the function x(t) is given by: 


2M 


em (t) = > aihi(t) (19) 
i=1 
The square of error norm for wavelet approximation in [1] is given by: 
ky 

— fy < —.—_ 2 

Therefore, 
1 
t) — eng (t)|| < O(— 21 
l(t) —ene(f)| < OC) (21) 


This means that error bound depends on level of resolution of Haar wavelets 
that is, error bound is inversely proportional to level of resolution of Haar 
wavelets. Therefore, when we increase the value of M, it yields the sure 
convergence of Haar wavelet approximation. 


6 Method of solution 


Consider the damped forced oscillatory equation (1). Assume that 


2M 
e'O =) ate). (22) 
i=1 
Integrating twice with respect to t from 0 to t, we get 
2M 
a’ (t) = 2"(0) + $7 a: Pi a(t), (23) 
i=1 
2M 
x(t) = 2(0)+ 5° a,P2,(t). (24) 
i=1 


Apply initial conditions and substitute the values of x(t), 2’(t) and x(t) in 
(1), we get, 

2M 

S © ai[Mhi(t) + rPi a(t) + SP2i(t)] = F(t) (25) 


i=l 
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where r, S, F and M are same as defined in Section 2. From here, wavelet 
coefficients a; are calculated and solution x(t) of equation (1) is obtained. 


7 Runge-Kutta method of fourth order 


Runge-Kutta method is famous numerical method for solving ordinary dif- 
ferential equations. Consider the second order ordinary differential equation 


aa) (26) 


By substituting dy = z, it can reduced to two first order simultaneous differ- 
ential equations 


d 
= 2 = Fle,y,2) (27) 
XL 
and 4 
ae 


) 
with initial conditions y(ao) = yo and z(ao) = 29. Starting at (0, yo, 20) 
and taking the step-sizes for x, y, z to be h, k,l respectively, the Runge-Kutta 
method gives, 


ky = hf (Xo, yo, 20); 29) 

l, = h(x0, yo, 20), 30) 
1 1 1 

kg =hf(to + 5 Yo 51 20 51), 31) 
1 1 1 

ly = ho(xo + 5s Yo 5 1, 20 51), 32) 
1 1 1 

k3 = hf (xo + 5 fs Yo 5 2 20 5!2)s 33) 
1 1 1 

lz = h@(xo + 5% Yo 5 2s 20 5!2)s 34) 

kg =hf (xo +h, yo + kz, 2% +13), 35) 

ly = hd(xo + h, yo + k3, 20 + Is). 36) 

Using above relations, we have 
1 
Yi = yor g (ki + 2ke + 2k3 + ka), (37) 


and 
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1 
Z1=zo+ g (a4 + 2ly + 2l3 + la), (38) 


To compute yo and zg, we simply replace x9, yo,2 by 2%1,4y1,21 in the 
above relations. Similarly by using above relations we compute 22, y2, 22, 
X35 Y35 Z3gercereceeee so on. 


8 Taylor-series method 


Consider equations (26), (27) and (28). If h be the step-size, y; = y(ao + h) 
and z1 = z(% +h). Then, Taylor’s algorithm for (26) and (27) gives 


/ h? N he wm 
yi = yo thyo + ayo + gp YO Tavs: (39) 
h? h3 
21 = 20 + hzo! t oy zo" 3! zol" seeeeeee (40) 
Differentiating (26) and (27) successively, we get y”, z’’, etc. So the values 
Yo, Yous Yous... and 20/, 207, 20", w+ are known. Substituting these values 


in above equations, we get y1, 21. Similarly, we have the algorithms 


h? h3 
yo =yi thy’ 4 ry yl 31 Yr HB ecesees (41) 
h? h? 
22 = 2% hz! t TI zy" 31 zy" FESO (42) 
Since y;, z, are known. we can calculate y;',y1", y1/", ...-- and 21) Zi ye ees 


Substituting these values in above equations, we get y2, 2. Proceeding in this 
way, we can calculate the other values of y and z step by step. 


9 Comparison of numerical solutions 


In this section, we compare the results of the present Haar wavelet method 
with two other numerical methods for the damped forced oscillatory problem. 
In order to verify the efficiency of Haar wavelet method in comparison to 
exact solution, Runge-kutta fourth order classical method and Taylor series 
method have been selected. For the Runge-kutta method, the step-size is 
1/32. For Taylor’s series method, step size is 1/32 and 7 terms are involved. 
Table-1 shows the numerical results from different numerical methods. Table- 
2 shows the errors arising from different numerical methods mentioned above. 
Further, graph in Figure 1 shows the comparison of graphical solution with 
the exact solution, obtained for J = 3 by (i) Haar wavelet method (ii) Runge- 
Kutta fourth order classical method and (iii) Taylor series method. 
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Table 1: Results from different numerical methods 


x(l)/32 | Exact solution | Haar wavelet | Runge-Kutta | Taylor series 
1 0.0004824193 | 0.0004707036 | 0.0004630874 | 0.0004671912 
3 0.0042356457 | 0.0042007228 | 0.0040510208 | 0.0040991982 
5 0.0114690029 | 0.0114111581 | 0.0109524412 | 0.0110920238 
7 0.0218950125 | 0.0218146369 | 0.0208783000 | 0.0211608344 
9 0.0352250351 | 0.0351226225 | 0.0335339697 | 0.0340203049 
11 0.0511707741 | 0.0510469174 | 0.0486211919 | 0.0493861034 
13 0.0694457517 | 0.0693011385 | 0.0658400073 | 0.0669763459 
15 0.0897667509 | 0.0896021596 | 0.0848906621 | 0.0865130140 


17 0.1118552180 


0.1116715126 


0.1054754806 


0.1077233303 


19 0.1354386196 


0.1352367450 


0.1273006949 


0.1303410858 


21 0.1602517482 


0.1600327242 


0.1500782254 


0.1541079134 


23 0.1860379714 


0.1858028875 


0.1732274013 


0.1787745029 


25 0.2125504200 


0.2123004292 


0.1973766160 


0.2041017526 


27 0.2395531092 


0.2392894218 


0.2213649076 


0.2298618522 


29 0.2668219898 


0.2665458670 


0.2452434581 


0.2558392936 


31 0.2941459241 


0.2938586718 


0.2687770046 


0.2818318053 


10 Conclusion 


Here, we used three numerical methods to approximate the solutions of 
damped forced oscillatory differential equation, and compared the results 
with exact solution. From above results, it is concluded that Haar wavelet 
method is simple, accurate and more efficient than other well known numer- 
ical methods for damped forced oscillatory differential equation. 
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Figure 1: Comparison of graphical solution of different numerical methods. 
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